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ABSTRACT 



Of 



Through the McDonald Observatory M Dwarf Planet Search, we have acquired 
nearly 3,000 high-resolution spectra of 93 late-type (K5-M5) stars over more than a 
decade using HET/HRS. This sample provides a unique opportunity to investigate the 
occurrence of long-term stellar activity cycles for low-mass stars. In this paper, we 
examine the stellar activity of our targets as reflected in the Ha feature. We have 
identified periodic signals for 6 stars, with periods ranging from days to more than 10 
years, and find long-term trends for 7 others. Stellar cycles with P > 1 year are present 
for at least 5% of our targets. Additionally, we present an analysis of the time-averaged 
activity levels of our sample, and search for correlations with other stellar properties. In 
particular, we find that more massive, earlier type (M0-M2) stars tend to be more active 
than later type dwarfs. Furthermore, high-metallicity stars tend to be more active at a 
given stellar mass. We also evaluate Ha variability as a tracer of activity- induced radial 
velocity (RV) variation. For the M dwarf GJ 1170, Ha variation reveals stellar activity 
patterns matching those seen in the RVs, mimicking the signal of a giant planet, and 
we find evidence that the previously identified stellar activity cycle of GJ 581 may be 
responsible for the recently retracted planet f (Vogt et al. 2012) in that system. In 
general, though, we find that Ha is not frequently correlated with RV at the precision 
(typically 6-7 m/s) of our measurements. 



1. Introduction 



The study of Ha activity for M dwarfs has provided a wealth of insight into the physics of 
low-mass stars. For young M stars with Ha emission lines (s o-called "active" M dwarfs), stellar 



activ ity has been shown to be tightly coupled with rotation (e.g. lBrowning et al.ll2010i : lReiners et al 



20121 ) . Furthermore, as these stars age and spin down, their Ha emission and variabili ty decrease 



significantly, offering a diagnostic for di stinguishing stellar populations of different ages (jWest et al 



2008 



2009 



Lopez-Santiago et al.ll2010l ). 



A bsorption line variability in inactive M dwar f s is often caused by pe r iodic stellar activity 



cycles (IKiirster et al. 120031; ICincunegui et ajj l2007al : IZechmeister et al.l |2009| ; iBuccino et al.l 12011 
Gomes da Silva et al.l |2012| ). Because such cycles may appear in radial velocity (RV) measure- 



ments (as evidenced by correlation between RV and activity tracers such as Ca II H and K emis- 
sion), potentially mimicking Keplerian planet signals (jQueloz et al.ll200ll ; llsaacson &; Fischerll2010l ; 
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Gomes da Silva et al.ll2012l ). they are of particular interest to planet search surveys. As a some- 
what happy coincidence, the long-term, multi-epoch spectral monitoring of stars conducted by RV 
surveys makes them uniquely sensitive to analogs of the 11-year solar cycle. These cycles may not 
appear in photometry if they manifest as variable heating of active chromospheric regions and not 
starspot modulation. 

Identification and characterizat ion of solar-type cycles will be instrumen tal in better under- 

) . Since ste llar cycles 



standing the solar magnetic dynamo ([Thompson et al.l l2003: W right et al 



are reflect ive of a star's mag netic field strength and variability (e.g. 



2011 



Brown et al 



structure (jOssendrij verl 120031 . and references therein), and long-term evolution (jWright 



2011 



internal 



20041 1. ex- 



ploring activity across all spectral types is an excellent way to understand how fundamental stellar 
properties vary, and how stellar magnetic fields are generated and maintained. 

In the Sun, the origin of the magnetic activity driving the observed 11-year cycle is believed 
to be the aQ dynamo maintained through differential rotation at the tachocline, the interface 
betwe en the radiative and convective layers of the Sun. The prevalence of similar cycles for FGK 
stars (jBaliunas et al.l Il995l ) confirms that such dynamos are common amongst solar-type stars. 
However, the t achocline lies at increasin g depths at later spectral types, disappearing altogether 
around ~ M4 (IChabrier &: Baraffdll997l ). Thus, a direct comparison between the stellar activity 
levels and cyclic behavior of solar-type and lower-mass stars is essential to characterizing the effect 
of stellar mass on internal architecture and the resulting magnetic activity. 

In this paper, we present the firs t systematic analysis o f stellar activity for the McDonald 



Observatory M Dwarf Planet Surv ey (|Endl et al 



Mount Wilson (jBaliunas et al.lll995l ) and HARPS (jLovis et al.ll201lf) have already conducted similar 



2003, 



2006). While the long-term surveys at 



studies for over 400 quiet FGK sta rs, analyses of M stars (e.g. lKiirster et al 



2003 



Cincunegui et al 



2007al : iGomes da Silva et al.ll2012l ) are significantly lacking, with less than 50 total stars examined. 
With 93 stars included in this paper, it is the largest such study to date for inactive M stars, and 
represents a substantial addition to the total collection surveyed. 

Here, we will use the flux in the Ha line as a tracer of stellar chromospheric activity. In 
areas of the stellar chromosphere where magnetic field lines influence the local convective behavior 
(so-called "active regions" or plages), Ha photons are emitted, resulting in an activity-dependent 
depth of the Ha absorption line. While this effect is more commonly observed in the calcium H 
and K lines, Ha is potentially a more suitable line for M stars because of their lack of flux near 
the calcium lines. Additionally, because Ha and Ca H and K are emitted from different depths 
in the chromosphere (and hence different distances from the tachocline), valuable information may 
be gained by comparing the results of M dwarf activity using both indices. Since Ha is a standard 
tracer of M dwarf activity, our measurements can easily be compared with previously published 
results. 

The paper is organized as follows. In Section 2, we describe how we acquired and reduced 
our data. We have identified a number of periodic Ha signals, which we discuss in Section 3. 
Additionally, we have computed the mean Ha flux levels of our targets, and include a detailed 
examination of how these average fluxes depend on stellar mass and metallicity. The broader 
implications of our results are discussed in Section 4, and we summarize our conclusions in Section 
5. 
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2. Observations and Data Analysis 



The McDonald Observatory Planet Search (e.g. Cochran h Hatzeslll993 ; Robertson et al.|[2012al fbh 
has monitore d nearby stars for R V variation induced by exoplanets with the Hobby-Eberly Tele- 
scope (HET. iRamsev et al.l Il998l ) since 2001. During this period, we have also surveyed 100 M 
dwarfs to explore the frequency of planets around low-mass stars. Our 100 targets are specifically 
selected to have low activity based on a lack of X-ray emission from the ROSAT All-Sky Survey 
( Hunsch et al.lll999l ). and should therefore represent the old, quiet population of nearby M stars. 



In order to better characterize our targets, we computed stellar masses using the lDelfosse et al 



(200 0j) K-band mass-luminosit y relation. Additionally, we calculate stellar metallicity using the 
Schlaufman fc Laughlinl (|2010l ) photometric metallicity calibration. We note that th i s calibr ation 
is largely consist e nt wi th earlier photometric metallicity estimates of iBonfils et aJj (|2005al ) and 
Johnson &: Appsl (|2009i ). but is less susceptible to systematically over- or underestimating stellar 
[M/H]. Where available, we have also obtained spectral subtypes from the Orion Spiral Arm Cat- 
alogue (jBobylev et al.1 120061 ) . In Figure CO we show the distributions of these stellar properties for 
our sample. 

After removing double-lined spectroscopic binaries (SB2s), the M dwarf survey has amassed 
2933 high-resolution spectra of 93 stars as of July 2012. All of the 93 targets have been observed 
at high cadence for a brief amount of time (typically ~ 5 observations over one week) to explore 
short-period RV variability, and are observed at least once per season to ensure sensitivity to long- 
period signals. Stars with potentially interesting RV signals have been observed several times each 
season. 



Our HET M dwarf spectra are taken with the High Resolution Spectrograph (HRS. iTull et al 

1998). Nearly all the spectra are taken at a resolving power R = 60,000, with a small number of 
bright targets for which we observed at R = 120, 000. We note that while our observations are 
taken with an I2 absorption cell in front of the slit for the purpose of obtaining precise RVs, no 
I2 absorption lines are present in the spectral region around Ha. Flat-fielding, bias subtraction, 
cosmic ray removal, and wavelength calibration are performed with standard IRAF0 routines, and 
we have divided our spectra by the blaze function to normalize the continuum. 

Although i t is common practice for RV surveys to trace stellar activity with the C a II H and K 
lines (ShK; e.g. iPaulson et al.1 12004 ; llsaacson &: Fischerll2010l ; iRobertson et al.ll2012al ). those lines 
are not accessible to HRS. Instead, we examine the Ha line, which appears i n absorption for our 
targets, but becomes increa singly filled in as chromospheric activity increases (jKiirster et al. I l2003l : 



Gomes da Silva et al. 



2011 



Figure [2] shows the Ha absorption line for GJ 270 at times of high 
and low emission. In addition to the issue of availability, analyzing Ha reveals stellar behavior 
which will not necessarily appear in S*hk since Ha activity a nd Ca H and K activity do not always 



corre l ate, as shown in previous co mparisons between tracers (jCincunegui et al.ll2007bl ; I Santos et al 



Gomes da Silva et al 



20101 
/ho- Following 



2011 



Kiirster et al 



We measure this chromospheric Ha "filling in" with an index 



(|2003l ). we define the index as 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 



-4- 



F l + F 2 



(1) 



As in iGomes da Silva et al.l (|201ll ) , we have taken i*H a to be the sum of the flux in a band of width 
1.6 A centered on the Ha A6562.828 A line. We find the line center with our fitted wavelength 
solution, making manual corrections for large stellar RV offsets where necessary. F\ and F 2 are the 
sum of the fluxes in the reference bands [—700 km s~ , —300 km s -1 1 (from the Ha lin e center) for 
band 1 and [600 km s _1 , 1000 km s _1 ] for band 2, as defined by iKiirster et al. I (j2003h . To ensure 
any variation observed is in fact from stellar chromospheric activity and not instrumental effects, 
we also measured a similar index for the Ca I A6572.795 A line, where -Fbai is the sum of the flux 
from [441.5 km s , 472.5 km s _1 ] around the Ca I line center, and using the same reference bands. 
Because the Ca I line should not be sensitive to stellar activity, the Ca I index serves as a control 
against which to verify Ha activity. 

For each stellar line index, we calculated an error bar using a method analogo us to the one em- 



ployed for the Mount Wilson Ca II Shk index we compute for our 2.7m spectra (e.g. lRobertson et al 



2012al ). Specifically, we multiplied the RMS in the continuum in the 0.5 A adjacent to the Ha or 
Ca I line by yN, where iV is the number of pixels in the line (i.e. <7f Hci = RM!3^ Ha±05 x y/N). 
Likewise, we multiplied the RMS scatter of each reference band by the square root of the number 
of pixels in the band, and added the errors in quadrature for the final line index. In equation form, 



+ 



'Pi 



+ a 



F 2 



F 1 +F 2 



(2) 



Rather than take the approach of IGomes da Silva et al.1 (|201ll ) and bin our observations, we assign 
these error bars so as to accurately assess the quality of each spectrum, and to best make use of 
our more sparsely-sampled data sets, which nonetheless show significant activity. 



3. Analysis 

We see time-dependent stellar activity for a number of stars from our RV survey. In addition 
to detailing this behavior, we have also examined the time-averaged overall activity levels of our 
targets, and related those levels to other fundamental stellar properties. We present these results 
in separate subsections below. 



3.1. Periodic Ha Activity 



We began our analysis by searching the time series In a measurements of each star for periodic 
activity. Since M dwarf stellar activity signals h ave been observed w i th periods as short a s weeks 
or months (as is typical of rotat i on per i ods, e.g. 



Kurster et a 



long as years (ICincunegui et al.l l2007al : IGomes da Silva et al 



2003: iForveille et al.l 12003 ) and as 
20121 ) . we searched a broad range 



of frequencies. We search fo r periodicity using the fully generalized Lomb-Scargle periodogram 
(jZechmeister &: Kiirsterll2009l ) due to its ability to handle uneven time sampling and individually- 
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weighted data points. In cases where our periodograms show significant power levels, we have 
estimated false alarm probabilities (FAPs) using three different me thods. The first two computa - 
tions, which we consider p reliminary estimates , use E quation 24 of [Z echmeis ter &; Kiirste r (2009) 
and the method outlined in lSturrock fc Scargld (|2010l ). For the lZechmeister fc Kiirsterl (|2009j) FAP 
estimate, we normalize the periodogram assuming gaussian noise, so our false alarm probability for 
a peak of power P is 

2P 



FAP = 1 - [1 - (1 



v- 



n) 2 ] 



M 



where N is the number of data points an d M = ^j- is an est i mate of the number of independent 
frequencies sampled. The FAP estimate of lSturrock fc Scargld (J201CJ), which compares the power of 
a periodogram peak to an expected distribution of power values based on Bayesian statistics, agrees 
well with the probabilistic calculation for our power spectra. Our final FAPs for p eriodic signals 



presen ted herein are calculated according to the bootstrap resampling technique of Kiirster et al 



(|1997l ). Briefly, the bootstrap FAP technique retains the time stamps of the original data set, then 
creates a series of fake data sets by assigning a random Ih q value (with replacement) from the set 
to each time stamp. The FAP is then taken as the percentage of resampled data sets which at any 
frequency give an equal or higher periodogram power value than the highest peak in the original 
data. 

Because stellar activity is neither strictly periodic nor perfectly sinusoidal, and will experience 
stochastic behavior (flares, etc.) regardless of any regular cycles, the FAP required to confirm 
a Ha activity cycle is somewhat higher than for a more rigidly perio dic phenomenon such as a 
Keplerian planet orbit. Previous studies (e.g. ICincunegui et al.ll2007al ) have considered FAPs as 
high as 0.5. While we see many signals in our sample with a FAP in the 0.5-0.1 range, we do not 
present them here, pending confirmation either through more dense time sampling (for short-period 
signals) or a second activity tracer such as the Na I D feature (for low-amplitude signals). Here, we 
only claim detection of a periodic signal if it has a bootstrap FAP at or near the 0.01 level, with the 
exception of GJ 581 (FAP ~ 0.1), for which we offer a confirmation of a pre viously observed cycl e 



(see below). This criterion is essentially equivalent to the one adopted by iBuccino et al.l (|2011 



The periodograms for our confirmed activity cycles are shown in Figure El and the individual FAPs 
are discussed in later sections. 

For each activity cycle, we have fit a sinusoid of the form Iua{P) = a o + a i sin(o;t + </>), where 
u) = 2p, P is the period of the cycle, and ao, ai, and <p are free parameters to set the In a zero 
point, amplitude, and phase, respectively. In Figured! we show the time-series Ha index of the six 
stars for which we observe periodic activity cycles and the fitted sine curves. Included with each 
plot is the time series of the Ca I index. We note that we have examined the periodograms of the 
Ca I index for each star in Figure HJ and do not see periodicity matching the signals seen in Ha. 
Furthermore, the level of variation in Ca I is considerably lower than we see in Ha, confirming that 
the cycles observed are in fact caused by stellar activity. 

While we show sinusoidal fits to the activity cycles shown in Figure (U it is clear from visual 
inspection that the cycles are not all perfectly sinusoidal in shape. Interestingly, the two most obvi- 
ously non-sinusoidal signals-GJ 552 and GJ 630 (see below)-both exhibit gradual declines and rapid 
rises in Ih q - This behavior is reminiscent of the Waldmeier effect seen in t he Sun, where higher- 
amplitude activity cycles are preceded by more rapid rises in activity (see, e.g. lCameron &; Schiissler 
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2008). While the details of the Waldmeier effect are the subject of active research, its qualitative 



properties can be explained by the variable nonlinear feedback of hydrodynamic turbulence and 



small-scale magnetic fields on the large-scale stellar magnetic field (jPipin Kosovichevll201ll ). If 
this is the case for our targets, we should expect to see varying shapes and durations of each cycle 
as we observe successive periods. 

We find that in cases of non-sinusoidal Iu a series, a better approach is to fit a Fourier series. 
In Figure O we show the results of fitting three terms of the Fourier series 

n 

Iu a (t) = a + ^2 a i sin(^it + 4>i) 

i=l 

for two of the cycles presented in Figure [H We set the initial guess for uj\ equal to the frequency 
identified in the periodogram of each star's Ihq index. For both stars, the best-fit frequency of each 
successive term falls very close to twice the frequency of the previous term, leading us to conclude 
that we are likely observing harmonics of the initially-observed activity cycle. The drawback to 
fitting Fourier series is that each new term adds three free parameters to the fit. For targets without 



a large number of observations (e.g. GJ 552, Figure 5(b)), we quickly lose the ability to constrain 



the parameters of higher-order terms, leading to over-fitting as shown in the Figure. 

Below, we discuss each observed activity signal in detail. Additionally, we have included the 
details of our fits to the data in Table [TJ 



3.1.1. GJ 270 and GJ 476 



With periods > 3 years, the Ih« signals of GJ 270, GJ 476, and GJ 581 (discussed below) could 
be considered the "stellar cycles" in our survey. Here, we see the advantage of utilizing the data set 
of an RV survey; the same observational strategy intended to reveal long-period giant planets gives 
us excellent coverage of very slow activity cycles, the durations of which have never been previously 
observed for M dwarf stars. 

At M* = 0.68MfD, GJ 27 (M2) is among the more massive M dwarfs in our survey. As first 



reported in lEndl et al.l ([2006), we see evidence for a long-period binary companion in the RVs of 
GJ 270, but at present its period is longer than our observational time baseline. The Ihq time 



series, shown in Figure 4(a) shows clear cyclic variation, with a distinct peak in the corresponding 
periodogram (Figure 3(a) ) at 2553 days. Our bootstrap analysis generated no false alarms in 10 4 
trials, thus giving an upper limit on the FAP of 10 -4 . We fit a sinusoidal model to this signal with 
a period of 2687 days and amplitude 0.00158 in Ihq; the longest-period cycle for which we have 



observed a full period. We include our fit to the data in Figure 4(a) 



Bobylev et al.l ([2006) list GJ 476 as subtype M4, although with a mass estimate of O.47M0, it 



would be among the more massive M4 stars. In Figure |4(b)| we show the 7h« time series, which 
again shows long-period variation. The power spectrum for GJ 476 peaks at 993 days, with a 
bootstrap FAP estimate again giving an upper limit of 10 -4 . Our fitted sinusiod gives a final 
period of 1066 days. The cycle has an Ih q amplitude of 0.00193, the largest of the periodic signals 
presented herein. 
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3.1.2. GJ 581 



Bonfils et al. 


2005b: 


Mavor et al. 


2009; 


Voet et al. 


2012) 



periodicity confirms the discovery of iGomes da Silva et al.l (|2012l ). It is important to note that our 



periodogram analysis ( Figure 3(c) ) shows a p e ak at 448 days with higher power than the longer- 
period signal claimed in lGomes da Silva et al.1 (120121 ) . However, since the longer-period peak leads 
to a marginally better fit-yielding an RMS of 0.00118 in 7ua versus 0.00123-we adopt the 1633-day 
peak as the true period. The 448-day peak disappears in the residuals to a 1633-day fit, indicating 
it is likely an alias of the fitted signal. 

As seen in Figure [3j the 1633-day signal of GJ 581 appears at somewhat lower power in 
our periodogram analysis. In our bootstrap FAP analysis, our randomly-resampled /jja indices 
produced higher power spectrum values in 1196 of 10 trials, for a FAP of 0.12. However, since 
our observed period is similar to the 1407-day period derived from the HARPS data, and because 
our observed /h» max i mum in December 2007/ January 2008 matches the Ha maximum found by 
Gomes da Silva et al.1 (|201lh . we consider our detection valid despite having a FAP of ~ 0.1. It is 



especially remarkable that we observe the cycle in Ha, while IGomes da Silva et al.l (J2012J) use the 
Na I D feature. Furthermore, the reversal of the minima/maxima between our Ha index and the 
HARPS Na I index confirms the anticorrelation between Ha and Na I for GJ 581 observed in the 



HARPS survey (jGomes da Silva et al. 



2011 



3.1.3. GJ 708 



With a mass M* = 0.77M Q , GJ 708 is at the upper limit of where the iDelfosse et al.l (120001 ) 
mass calibration is valid. While the Orion Spiral Arm Catalog lists it as spectral type Ml, its mass 
indicates GJ 708 is most likely a late K dwarf. Our RV series shows a monotonic decreasing trend 
indicative of a binary companion, for which we are currently unable to estimate a period. In many 
ways, then, GJ 708 is quite similar to GJ 270. 

GJ 708 is unique among the objects presented here because it exhibits a long-term trend in 
Ha along with a periodic signal. We note that we have fitted and removed a linear slope from the 
/na series prior to creating the phase plot shown in Figure [4(d)] Once removing the trend, the 
periodogram on the residual Ih q values (Figure [3(d)[ ) shows significant power at 296 days, although 
we note that the peak appears regardless of whether the trend has been removed. In 10 4 bootstrap 
resamplings of the residual /jia values for GJ 708, we observed 282 false alarms, for a FAP of 0.03 
on the 296-day period. It is also important to point out that, while the period observed for GJ 708 
is close to the 1 year alias, neither the periodogram of our time sampling (the window function) 
nor the Ca I periodogram showed periodicity near 296 days, eliminating the possibility that the 
observed signal is caused by our sampling. 

The 296 day period of GJ 708's activity cycle is too long to be the rotation period of a typical 
old dwarf star (normally weeks or months), but is also somewhat short to be reminiscent of a long, 
solar-type cycle. On the oth er hand, it is quite simil ar to the ~ 442 day period observed for the 



M dwarf Proxima Centauri (jCincunegui et al.l l2007al ) . It is possible, then, that these two stars 
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represent the first examples of a new class of intermediate-duration activity cycle that occurs in 
low-mass stars. Alternatively, it may be the case that the linear trend we observe is the star's true 
"activity cycle," and that the 296-day signal is an intermediate-peri od cycle that exists a longside the 
longer signal. Such sub-cycles have been observed for the Sun (e.g. 
similar physics may be at work here. 



Tan fc Chengj|2012l ). suggesting 



3.1.4. GJ 730 



The periodogram for the time series /ho of GJ 730 (Figure 3(e) ) shows interesting peaks at 



3.06 days and 993 days. While the longer periods of the other Ih q signals in this paper make the 
longer period peak seem intuitively more likely to be the true signal, the 3-day signal continues to 
increase in power as we acquire additional data, while the 993-day peak remains constant. While 
many of our stars exhibit short-term variability, only GJ 730 displays a coherent signal of such 
statistical significance; our bootstrap FAP estimate gives an upper limit of 10 for the 3-day peak, 
making it much more significant than the long-period peak. We note that the 993-day peak is no 
longer present in the residuals around the 3-day fit, and also that the 3-day peak disappears from 
the residuals to a 993-day fit, giving further evidence to the two signals being aliases. Because it 
is possible for an alias to display stronger periodogram power than a true signal, and because the 
longer period is much more typical of stellar activity cycles, we will not make a definitive argument 
as to which is the true period. We include fits to both periods in Table [Vj and show phase plots of 



both signals in Figure 4(e 



If the true period for GJ 730 is 3 days, it is by far the shortest periodicity presented herein. 
While 3 days could concievably be a rotation period for a very young, active star, it is almost 
certainly too short for a main-sequence M star, although we are unaware of any v s in i measur ements 
which might confirm or reject the rotation hypothesis. Hipparcos photometry (|ESAlll997l ) is not 
sampled densely enough to be sensitive to a 3-day period. At spectral type M2, with a mass of 
O.62M0 and [M/H] = 0.014, GJ 730 appears more or less typical of our sample, leaving no indication 
as to why such an unusually short periodicity should be present. Its average Iu a value of 0.064 is 
also close to the overall average for our targets, so the star is not especially active. 

We note that 3 days is common for the periods of hot Jupiter planets. Previous results 

L) have shown evidence that hot Jupiters may induce 



(Shkolnik et al 



2008; 



Pillitteri et al. 



2011 



star-planet interaction (SPI), with measurable periodicity in stellar activity tracers matching the 
period of the planet's orbit. Our RVs do not currently show evidence for planetary companions 
to GJ 730, and the RMS scatter of 7.4 m/s shows definitively that any hot Jupiter planet would 
need to be in an orbit which is highly inclined relative to the line of sight. While the inclination 



requirement in combinati on with the inherent paucity of hot Jupiters around M stars (jEndl et al 



2006 



Johnson et al.ll2012l ) makes SPI an unlikely explanation for the 3-day activity signal observed 



for GJ 730, we present it as an interesting possibility. 

On the other hand, if the ~ 1000 day signal is the true period, GJ 730 fits nicely with GJ 270, 
GJ 476, and GJ 581 as a star displaying solar-type activity cycles. Further study will be required 
to determine the nature of this star's periodic activity. 
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3.1.5. GJ 552 

GJ 552 is an extreme example of non-sinusoidal behavior in Ha activity cycles, as evidenced by 
the sharp upturn in Jho observed between July 2006 and April 2008. As mentioned previously, while 
we recognize the inadequacy of a sine wave fit, we do not have a sufficient number of observations 
to fit a long Fourier series. Furthermore, we see in Figure |5(b)| that our different fits offer quite 
dissimilar estimates for the period of the cycle. While it is possible we have recorded close to a 
full period, some models result in a period far exceeding our observational time baseline. For this 
reason, while we believe the behavior seen in our Ha index is likely periodic, we do not claim a 
specific period, and do not include a model fit in Figure |4(f)[ 



3.2. Long-Term Activity Trends 

For a second subset of our target stars, we see long-term slopes or curvature in Iu a , but are 
unable to confirm any periodicity for these objects because our observational time baseline is shorter 
than the duration of any periodic signals. We show the time-series Iua data for these stars-GJ 16, 
GJ 521, GJ 96, GJ 3023, GJ 3801, GJ 611.3, and GJ 630-in Figured along with the corresponding 
Ca I index. The details and statistics of our fitting are included in Table Q] as well. We separate 
the stars in Figure [6] into two groups: stars showing monotonic slopes, and objects for which we 
observe some curvature or turnaround. 



3.2.1. GJ 16, GJ 521, GJ 96, and GJ 3023 



GJ 16, GJ 521, GJ 96, and GJ 3023 all display a linear trend in J H a over our entire observational 
time baseline. To ensure the statistical significance of these slopes, we perform a linear least-squares 
fit to each time series, and calculate the Pearson correlation coefficient. We include these fits in 
Figures |6(a)||6(d)[ In all four cases, the correlation coefficient indicates the relations are highly 
significant (with the probability of a null result p « 0.01). 



3.2.2. GJ 3801, GJ 611.3, and GJ 630 



GJ 3801, GJ 611.3, and GJ 630 all show significant variability, but are poorly fit by linear 
regression. For GJ 3801 and GJ 611.3, we have compared the quality of linear and quadratic fits 
with an F-test using 



poly 



(N 



3)(X 2 



slope 



^poly)/ X 



2 

poly 



as described in IGomes da Silva et al.1 (|201ll ). Here, Xsiope * s ^ ne chi-squared value of a straight- 
line fit, while Xpoiy * s the chi-squared of a quadratic fit. We find P(F po i y ) = 0.04 for GJ 3801 and 
F(F po i y ) = 0.01 for GJ 611.3, indicating the quadratic fits are a statistically significant improvement 
over straight-line models. Our fits are included in Figures 6(e)||6(f) In GJ 630, the recent dramatic 
upturn in Ihq looks very similar to that of GJ 552. If the Ha activity seen for GJ 630 does prove 
to be another highly non-sinusoidal periodic signal, its period will easily exceed 10 years, making 
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it a direct analog to the 11-year solar cycle. Since the behavior of In a is approximately linear prior 
to this rapid upturn, we have not included a fit for GJ 630, pending future observations to more 
accurately characterize the shape of the curve. 



3.3. RV correlation 

Among our motivations for examining the stellar activity of our targets was to evaluate to 
what level, if any, chromospheric variability influenced our RV measurements. To this end, we have 
searched for correlations between RV and 7ho f° r our entire sample. In cases where we see evidence 
for substellar companions (planets or brown dwarfs) in the RVs, we have subtracted those signals 
and examined the residual velocities. 

Over all, we see little corre l ation between Ih« and RV. We find this result unsurprising for 



our data; iGomes da Silva et al.l (120121 ) report correlations between stellar activity indicators and 
RV typically occur at or below the 5 m/s level. Such an effect will doubtless be significant for 
our sample following the HET/HRS upgrade, which will yield a throughput gain of approximately 
2 magnitudes at our R = 60,000 setting, dramatically improving the RV precision for our faint 
M dwarf targets. However, with our current RV precision limited to ~ 6 m/s for the majority 
of our sampl e, it is unde r stand able that we see relatively few correlations between /jja and RV. 



Furthermore, lLovis et al.l (|201ll ) observe a decrease in correlation between RV and -Rjjk w ith de- 
creasing T e fj, so it is likely we will see relatively fewer M stars with RV-Ha correlation regardless 
of precision. 

On the other hand, we do see two examples-GJ 1170 and GJ 3801-for which RV is anticorre- 
lated with Jhq at a level higher than 5 m/s. We show our measured RVs and the corresponding 
Ha indices for these stars in Figure [71 For GJ 1170, we derive a Pearson correlation coefficient 
of -0.48 for the relation between RV and Ift a , yielding a probability p = 0.005 of no correlation. 
Similarly, for GJ 3801 the Pearson correlation coefficient is -0.51, indicating the probability of no 
correlation is just 0.03. 

GJ 1170 is a particularly interesting the Ha emission is reflected in the RVs, resulting 

in a signal which initially appeared to be a Keplerian planet orbit with RV amplitude K ~ 20 m/s. 
Assuming a stellar mass of M* = 0.52M Q , this signal corresponds to a planetary minimum mass 
of approximately 0.9 Jupiter masses. We therefore conclude that, while In a does not frequently 
correlate with RV for M dwarfs at a precision level higher than 5 m/s, there does exist a small 
subset of stars where Ha activity corresponds to large RV shifts. With such a small number of 
objects, though, we must concede the possibility that these stars actually have giant planets causing 
the observed RVs, and coincidentally exhibit activity cycles with periods similar to those of the 
planets (as with Jupiter and the solar cycle). 



3.4. Stellar Activity Levels 



In addition to examining time-series activity of individual stars, we have also considered the 
mean activity levels of our targets relative to the overall sample. To do so, we average our Ha index 



- 11 - 



over all observations of a star. 

Based on their observed positive correlation between time-averaged Ha index and V — I color, 
Gomes da Silva et al.l (|201ll ) report that mean activity is a function of temperature. We attempted 
to confirm this correlation with our sample, but find considerably larger scatter than previously 
observed, as shown in Fig ure El While the relation hold s relatively well within the 1.6 < V — I < 2.8 
color range covered by the lGomes da Silva et al.l (|201ll ) sample, it is much less suitable outside that 
regime. While we looked at mean I-R a as a function of several color indices, we found that no color 
was a satisfactory tracer of Ha. 

We considered the possibility that stellar mass might be a better predictor of mean Ha activity 



than temperature or color. In Figure 9(a) , we show our time-averaged Ha index as a function of 



stellar mass. The resulting anticorrelation is obvious, and the scatter is greatly reduced from that 
of the color relations. It is apparent that Ha emission is much more tightly connected to mass than 
temperature. Performing a linear least squares fit to the data, we find the relation 



(i- H ) = 0.094 - 0.050^ 

v ' M 

Our fit gives la uncertainties of 0.003 on the slope and 0.001 on the intercept. 



(3) 



We explored the residuals to the fit above to determine whether additional dependencies exist. 
As before, we see no convincing correlation to any color index. On the other hand, stellar metallicity 
does appear to play a role. Figure |9(b)| shows the residual Ha index after subtracting the Ha -M* 
relation. The resulting correlation is a clear indicator that stellar [M/H] is the primary contributor 
to the scatter around the mass-activity fit. Again performing a linear least-squares fit to our data, 
we find 



(^Ha)r 



7.6 x 10~ 4 + 9.7 x 10~ 3 [M/H] 



(4) 



with la uncertainties of 10 3 on the slope and 3 x 10 4 on the intercept. 

Before drawing conclusions as to the physical interpretation behind Figure [9l it is important 
to understand how the stellar continuum changes as a function of mass and metallicity. Since 
In a is normalized by the continuum adjacent to the Ha line, it is possible that the observed trends 
are due to changes in continuum flux rather than Ha emission. We considered this possibility by 
converting our Ha index to an equiva lent width (EW) using a transformation basically identical 
to the one outlined in Appendix A of IZechmeister et al.l (120091 ). then to the quantity [Ln a / L^oi], 
the logarithm of the luminosity in Ha, normalized by the star 's bolometric luminosity. For this 
transformation, we used the method of IWalkowicz et al.l (|2004l ). which calibrates the ratio of the 
continuum near Ha to the bolometric luminosity as a function of color for M dwarf stars. Briefly, 
we use the transformations 



EW Ha = AA Ha (l - ^/h«) 



(5) 



[Ln a /L hol ] = xEW Ha 



(6) 
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Where Ac = 8.75 A is the width of each reference band used to compute In a , Ah q = 6562.808 
A, AAho = 1.6 A is the width of the spectral wi ndow in which we sum t he Ha flux, and x is the 
luminosity scaling factor, calculated according to lWalkowicz et al.l (J2004J). In Figure [TUl we again 
show Ha emission as a function of stellar mass and metallicity, but measure Ha emission with 
[L Ha /L ho i] instead of I Ha . 

Perhaps surprisingly, both of the trends seen in Figure [9] are reversed in Figure [TUl Our fit to 
[Lua/Lboi] as a function of stellar mass gives 



M* 

([WA>oi]> = -4.17 + 0.87-4 

M 



(7) 



with residuals anticorrelated with stellar metallicity. A linear least squares fit to Figure 10(b) 
yields 



res 



-0.019 -0.27[M/H] 



(8) 



Our la uncertainties on the slopes and intercepts, respectively, are 0.08 and 0.04 for Equation [71 
and 0.01 and 0.04 for Equation [HJ 

At first, the results of Equations [3M] and [7J[8] appear to present conflicting information regarding 
the dependence of stellar activity on mass and metallicity. However, we suggest a conceptually 
simple hypothesis consistent with all four relations. It has bee n firmly established that the stellar 



continuum near Ha var ies along with bol ometric luminosit y ( Hal 



1996 



Walkowicz et al. 2004 



Cincunegui et al.ll2007bl ). According to the lCincunegui et al.l (l2007bl ) -Fa6605 versus B — V relation, 
we expect the continuum around Ha to differ by up to a factor of 5 across the range of B — V typical 
of our sample. Thus, since [L^ a / L^ i] spans a factor of ~ 4 over our sampled range of stellar mass, 
we would expect a ~ 80 percent decrease in In a over the same range (assuming AIu a ~ A ^° 5 ). 
The observed In a decrease of ~ 60 percent is fully consistent with this expectation, considering 
the significant uncertainty in the Cincunegui relation at high B — V. We therefore conclude that 
chromospheric Ha emission is increasing with stellar mass, but the effect does not appear in our 
original Ha index because of the increasing continuum flux. 

Understanding the metallicity dependence for Ha emission requires a similar deconstruction 
of stellar luminosity/continuum effects. When we subtract a linear fit to the 7ho - M* relation, 
the effects of the varying stellar continuum should be removed, suggesting higher-metallicity stars 
are more active. However, verifying this demands explaining the decrease in [Ln a /L^ \] with 
[M/H]. Previous observations of M dwarf activity ha ye shown that more a ctive stars are more 



luminous, possibly due to an increase in stellar radius (jBochanski et al.ll201ll ). Furthermore, stars 
with high [M/H] will be more luminous at a given B — V or V — K color than metal-poor stars 



due to line blanketing of blue light in the more metal-rich objects. Figures [9(b)1 and [10(b) 



arc 



therefore consistent with the argument that metal-rich M dwarfs display higher Ha emission, but 
their increased activity and opacity in the blue cause L^ \ to increase faster than Lu a , causing the 



negative slope seen in Figure 10(b) 



Earlier studies of M dwarf activity (e.g. IWest et al 



2008 



Browning et al.l |2010| ) emphasize 



the importance of considering the various subtypes of M stars, particularly around M4, where 
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stars typically become fully convective. Figure [TT] shows the time-averaged [Lho/Ajoi] values for 
objects of known subtype in our sample. In the Figure, the error bars represent the RMS scatter 
within a given bin, rather than measurement error. The decrease in [Ln a /Lb \\ towards later types 
essentially reflects the observed mass dependence. However, we see more clearly in Figure [11] the 
tendency for more active bins to have more scatter as well. 



4. Discussion 



4.1. Stellar Cycles 



Our observed stellar activity cycles indicate that even for an old, quiet stellar population, 
magnetic activity may appear at essentially any period. The long-period subset of stars shown in 
Figure [U particularly GJ 270 (P ~ 7.5 years), demonstrate stellar periodicity on the order of the 
11-year solar cycle is far from unique; stars with cyclic activity of period > 1 year appear in at least 
5% of our sample . In future work we will include additional tracers (e.g. the Na I D feature, see 
Diaz et al.ll2007al ). both to increase sensitivity to additional activity cycles, and to fully evaluate 
the impact of stellar variability on our RV survey. The addition of other tracer indices will almost 
surely increase the number of candidate signals among our targets. 

While the period s of the stellar cycles ob served in FigureH]fall towards the low end of the period 
distribution found bv lBaliunas et al.l (| 19951 ) for solar-type stars, our results are entirely consistent 
with the Mount Wilson survey. While the Baliunas distribution peaks between ~ 5 — 12 years, 
there are a significant number of detections between 1 and 5 years. Furthermore, our observational 
time baseline, while extensive, is still too short to properly characterize stellar cycles with P > 10 
years. It is quite likely that the six stars shown in Figure [6] are examples of partial detections of 
such periods. Clearly, it is essential that both the durations and number of targets be increased for 
M dwarf activity surveys for a proper comparison to FGK stars. 



Lovis et al.l (|201ll ) find activity cycles with amplitude 0.04 or higher in -Rjj K for 61 ±8% of quiet 



FGK stars. At present, the perce ntage of M dwarfs with obse rved magnetic cycles is significantly 
lower. Using the Na I D feature, iGomes da Silva et al.l (|2012l ) find periodic signals for 5 of 27 M 
stars surveyed, for a rate of 19%. With 6 periodic signals found in 93 stars, our resulting "cyclic 
fraction" is even lower, at 6.5%. If we assume the 7 stars for which we observe long-term trends 
are exhibiting activity cycl es, our rate increases to 14%, which is more in agreement with the 
Gomes da Silva et al.l (|2012l ) value. It is possible that the difference in the periodic fractions is due 
to the relative sensitivites of the activity indices. Because the sodium feature is emitted deeper in 
the chromosphere than Ha, where the local density is higher, a change in the magnetic field should 
produce relatively more emission in sodium than hydrogen, making the amplitude of a given stellar 
cycle higher in the Na I index than in In a . Our overall sample size, and the number of observed 
stellar cycles and trends in Ha will allow us to make a thorough comparison between the two indices 
in future work, as we will have a robust set of comparison cases to establish whether the Na I D 
feature is simply more sensitive to weaker signals, or in fact probes different stellar physics. 



In any case, current data suggests M dwarfs are less likely to display measurable magnetic 
periodicity than earlier spectral types. A possible explanation for the discrepancy lies in the lack 
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of a tachoc line for many M stars. At a spectral subtype of M4 and below, stars become fully 
convective (jChabrier Baraffdll997l ) and therefore do not possess a solar-type dynamo, which is 
driven by differential rotation between the radiative and convective layers. If the solar cycle and 
its analogs are driven by this dynamo, then we should expect to see only low-amplitude cycles or 
quasi-periodic oscillations for fully convective stars. Indeed, of the stars we identify a s having cycles 



or tren ds, only GJ 476 and GJ 581 have spectral subtypes later than M3 according to lBobylev et al 



(|2006l ). Of those two, GJ 476 has an estimated mass of 0.47MQ, which is massive enough to have 
a tachocline despite its listed type of M4, and the activity cycle of GJ 581 may be connected to its 
planetary system (see below). Furthermore, Figure [TT] shows a distinct dropoff in [L-n a /L^ \] around 
spectral subtype M3-M4, suggesting the fully convective stars in our sample are characteristically 
less active. If continued study shows stellar activity cycles in M stars to be confined mainly to 
subtypes M0-M3, it will be an intriguing piece of evidence in favor of the tachocline dynamo as the 
origin of stellar magnetic cycles. 



The period of the In a cycle seen for GJ 270 is nearly equal to that of GJ 752 A (IBuccino et al 



201 ll ). making it one of the two longest M dwarf activity cycle s observed to date. It is worth 



notin g that both of these stars are members of binary systems (|Endl et al.l 120061 ; IBuccino et al 



201 ll ) . Among the stars in Figures 0] and [61 we see evidence in our RV s for a brown-dwarf b inary 
companion to GJ 708, and GJ 581 hosts a known multi-planet system. iMorgan et al.l (|2O120 have 
shown that M dwarfs with close binary companions display increased magnetic activi ty, and the 



SB2GJ 375 has b een shown to exhibi t a lon g-period activity cycle in the Ca K index (jDiaz et al 



2007b! ) . Likewise, iKrejcova Budajl (J2012J) see tentative evidence of increasing Ca II emission 
for stars in known exoplanet systems. Interestingly, of the binary /planetary systems for which 
stellar cycles have been observed, none of the stellar activity periods match the companion orbital 
periods or their harmonics. Further study will be required to definitively establish or reject a 
physical connection between more distant/less massive substellar companions and stellar activity 
cycles, but it is worth noting that a number of our candidates exhibit both, as do those of previous 
studies. 



4.2. Magnetic Activity and RV Planet Surveys 



Aside from the obvious value to the study of low-mass stellar physics, the understanding of M 
dwarf activity cycles has important implications for RV planet surveys. It has been well established 
that periodic RV signals can be caused by stellar chromqspheric variability, as traced by photometry 



or acti vity indicators such as Ca H and K (IQueloz et al.ll200ll : lForveille et al.ll2009l : lRobertson et al 



2012al ). While activity- induced RV signals are typically expected to have amplitudes of just a few 
m/s, Figure [7] shows that exceptionally large-amplitude effects will occasionally appear. In general, 
though, these results will be increasingly important as M dwarf RV surveys approach precisions 
of 1 m/s or less. GJ 581 stands as an excellent example. We notice a striking similarity between 



our periodogram of the GJ 581 Ha in d ex (F igure 3(c) ) and the periodogram of the residual RVs 
to a 4-planet fit shown in IVogt et al.l (|2010l . Figure 3, panel 5), particularly at the ~ 450- and 
~ 1600-day peaks. The 433-day period of the unconfirmed planet f is very close to our observed 
alias at 448 days, as is the candidate 400-day planet identified in iGregoryl (|201ll ). While fitting 
a Keplerian orbit to an activity-induced RV signal would essentially serve to correct the RVs for 
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stellar activity, if the true period of GJ 581's In a cycle is 1300-1600 days, then removing a 433-day 
orbit would lead to incorrect residuals, complicating the analysis of any additional planets. Clearly, 
any RV detection of low-mass planets around M dwarf stars should be accompanied by a careful 
analysis of the corresponding stellar activity indices. 



4.3. Mean Activity Levels 



Our results for the comparative activity levels between the stars in our sample are interesting 
in the context of previous surveys of more active M stars. Figure [TT1 s hows essentially identical 
behavior to the findings o f IWest et al.l (120041. F igure 5) and Kruse et al.l (|2010l . Figure 7). We note 
a slight discrepancy from iMorgan et al.l (|2012l . Figure 12), though, who do not see a decrease in 
[Lna/^boi] until spectral type M4. Whether this behav ior extends to solar-type stars is some- 
what unclear based on current results. lLovis et al.l (|201ll ) report that most G stars have log -Rjjk 
around —5, while stars of K type cluster around logi?^ K ~ —4.7, suggesting an inversion of the 
relation observed for M stars. These results could still be reconciled, since stellar bolometric lu- 
minosity increases at earlier spectral types, so some normalization may be required to properly 
understand how mean activity levels vary across the entire main sequence. Of course, it is known 
that Ha behaves differently f rom Ca H and K, especially for cool dwarfs such as those in our s ample 
(e.g. lCincunegui et al.ll2007bl ). A more harmonious result can be seen in lBaliunas et al.l (|1995l ). who 
plot mean log S*hk against (B — V) color for FGK stars, finding an increasing trend (and scatter) 
similar to Figure 

The data presented here are perhaps more compelling when shown as a function of stellar 
mass, though. Because our targets are pre-selected to b e "quiet," and do not show Ha in emission, 
they should all be older than their "activity lifetime" ([West et al.l 120081 ) . Figure 10(a) seems to 



suggest that for these quiet stars, the dynamo responsible for driving the residual magnetic actvity 
following the spin-down phase is more effective in more massive stars, which presumably are not 
fully convective. Again, it is possible that the presence of a more significant radiative zone in the 
interiors of the more massive stars leads to the generation of a solar-type magnetic dynamo at 
the interface of the convective and radiative layers, rather than the rotation-dominated magnetic 
activity of younger, highly active M dwarfs. 

As explained above, our data indicate that metal-rich stars tend to be more active than average 
at a given stellar mass. One possible cause for this relation is that our measured [M/H] is-in a 
statistical sense-a tracer of age, and we are simply ob serving the previously-confirmed conne ction 
between stellar age and chromospheric variability (e.g. lWest et al. 2008; Bochansk i et alj|2011f). T o 



test this hypothesis, we have obtained space velocities for our sample from lBobvlev et al.l ((2006). 
In Figure [T2l we plot the mean galactic velocity dispersions of various metallicity bins in the UVW 
coordinates, excluding very high and low metallicities due to small numbers of stars. In general, 
we see that low-metallicity stars have higher velocity dispersions in all three coordinates, lending 
support to the proposal that the metal-rich stars in our samp le are younger. Further evidence of the 
age-metallicity relationship is discussed in lwest et al.l ( 2008 ). The fact that we see higher average 
/Ha values towards increasing [M/H] is therefore likely a reflection of the younger members of our 
sample being more active. Again, all of the stars discussed herein are old, and should have evolved 
past their early active stages. It is remarkable, then, that the effects of stellar age are still visible in 
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the activity levels of these objects. It is worth pointing out that, w hile the UVW veloc ities suggest 
a spread in age among our objects, the probability distributions of iReddv et al.l (|2006l ) indicate all 
of our targets are members of the Galactic thin disk with > 99% probability, so our observed mass 
and metallicity/age dependencies may or may not apply as cleanly to stars in the thick disk or 
halo. 

The implications of this study are potentially encouraging for the prospects of current and fu- 
ture near-infrared (NIR) RV surveys for planets around late M stars. The old, low-mass stars which 
will make up the majority of the targets should have relatively low chromospheric activity levels, 
and will be less likely to exhibit activity cycles. Furthermore, we have shown that the stellar cycles 
which may be confused with planetary signals appear in Ha, which is much mor e likely to be acces- 



sible to a NIR spectrograph than Ca H and K. The CARMENES spectrograph (jQuirrenbach et al 



201 ll ). in particular, will be able to acquire Ha. Therefore, problematic signals may be properly 
identified and corrected. While the Habitable Zone Planet Finder on HET (jMahadevan et al.l l2010) 
will not have wavelength coverage at Ha, our results illustrate the necessity of identifying a suitable 
activity tracer in the NIR for RV surveys. 



5. Conclusion 

We have measured and analyzed time-series Ha emission for ~ 11 years of spectra for 93 
low-mass stars. Our primary results may be summarized as follows: 

1. We find strong periodic Ha variability for 6 stars in our sample, and long-term trends 
or curva ture for 7 others. Among th ese, we confirm the discovery of the activity cycle for GJ 581 



found bv lGomes da Silva et al.l (|2012j ). Our observed stellar cycles are analogous to those previously 



observed for FGK stars. 

2. While activity-induced RV signals are relatively rare in our sample at its current precision 
level, we identify two stars-GJ 1170 and GJ 3801-for which our Ha index is anticorrelated with RV, 
indicating activity cycles causing RV variations of nearly 20 m/s. Additionally, we see evidence that 
the signal of the purported planet f in the GJ 581 system may in fact be caused by an alias of the 
~ 1600-day stellar activity cycle, further emphasizing the need to consider stellar RV modulation 
at the < 3 m/s level of precision. 

3. The mean Ha activity levels of our targets correlate with stellar mass. This result is 
qualitatively consistent with previous observations of both solar-type and low-mass stars, although 
the correlation has more often been expressed in terms of stellar temperature or color, rather than 
mass. We confirm that [Ln a /L^ \] decreases towards later stellar subtypes, a phenomenon that 
appears to hold true for M stars regardless of whether or not they are in the "active" phase. Lower 
[LHa/^boi] values and an apparent lack of observable stellar cycles for fully convective stars may 
be related to the absence of a tachocline for those stars. 

4- In addition to the mass- activity relation, we see a metallicity dependence on average 
Ha emission. Lower average galacitc velocity dispersions indicate the more metal-rich stars in 
our survey are in fact younger. Even though all the stars in our sample have aged beyond their 
"active" phase, it appears some residual age effect still contributes to their mean activity levels, 
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and that the younger of our targets are still more active at a given stellar mass. 
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Fig. 1. — The McDonald Observatory M Dwarf Planet Search is a dedicated long-term survey of late-type stars in the Solar neighborhood. 
Here, we give histograms of the stellar masses, metallicities, and spectral subtypes of our targets. 
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Fig. 2. — For quiet M stars, the Ha line traces stellar magnetic activity through the variable 
"filling in" of the absorption line with photons emitted from magnetically active regions of the 
chromosphere. Here, we show the Ha absorption line of GJ 270, in states of high (red) and low 
(blue) chromospheric emission. The dashed vertical lines indicate the 1.6 A window in which the 
Ha flux is calculated. 
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Fig. 3. — Fully generalized Lomb-Scargle periodograms for the six stellar activity cycles detected in our M dwarf sample. Below each 
periodogram is the power spectrum of our time sampling (the window function). The horizontal lines represent the power levels corresponding 
to a false alarm proba b ility (FAP) of 0.5 (dotted line), 0.1 (dashed line), and 0.01(dash-dotted line), as calculated from Equation 24 of 
Zechmeister &: Kiirsterl ([2009). We note that these FAP levels represent a preliminary estimate, and our formal FAP values are obtained 
through a bootstrap analysis, which we give in Section 3. 
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Fig. 4. — We have searched the time-series Ha variability of our entire M dwarf survey for periodic activity, particularly analogs to the solar 
cycle. The stars in this Figure show definitive evidence of periodic activity in the Ha line. The data have been phase-folded for shorter 
periods for clarity. The red curves show our best fit to the data. For GJ 708, we have also removed a linear trend. In the case of GJ 730, we 
include two possible fits (see Section 3.1.4). We do not include a fit for GJ 552 because we are unable to satisfactorily constrain the period 
of the observed cycle. Below each Ha series is the corresponding Ca I index, which we use as a control. 
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Fig. 5. — Several of our observed activity cycles exhibit non-sinusoidal 7h« variability. Here, we 
show the results of fitting Ha activity cycles with Fourier series instead of sine curves for a GJ 
270 and b GJ 552. We see that our fits improve with the inclusion of higher-frequency harmonics, 
but that for few data points (as seen in b), we cannot constrain enough free parameters to justify 
adding a large number of terms. 





(g) 

Fig. 6. — In addition to clearly periodic signals, we observe seven stars with long-term linear trends or curvature in In a • In this Figure, we 
show the time-series Ha indices for those stars. Below each Ha series is the corresponding Ca I index, which we use as a control. In a-d, we 
show the linear least squares fit to the data in red. For e and /, the red curve indicates our best-fit quadratic function. 
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Fig. 7. — Previous results suggest the effects of stellar activity on RV measurements should be at 
or below the 5 m/s level, but we find two examples of activity- RV correlation at larger amplitudes. 
Here, we show radial velocity measurements (Top) and their corresponding Ha indices (Bottom) 
for GJ 1170 (a) and GJ 3801 (b). The red curve over the GJ 1170 RVs shows a Keplerian fit 
to the data. The RV variation appears to be caused by a stellar activity cycle, as seen by the 
anticorrelation between RV and /ho- 
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Fig. 8. — If stellar activity varies as a function of T e ff, the average Ha index should correlate with 
V—I color. In a, we show our time-averaged Ina as a function of V —I for our entire sample, b shows 
the sa me data, but restricted to the same range in V — I as covered by the iGomes da Silva et al 
(l201lh sample. 
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Fig. 9. — Since mean Ha activity does not correlate well with color (and therefore presumably 
T e ff, we investigate whether there are other stellar parameters which better predict Ha behavior. 
a. shows our time-averaged Ha index as a function of stellar mass for our sample. The red line 
shows our linear least squares fit to the data. In b, we give the residuals to the fit shown in a, 
plotted as a function of stellar metallicity. The red line shows our linear least-squares fit to the 
data. 



- 30 - 



; i 1 i 1 i 1 r 

-3.2 - • 




-0.4 - 

J i I i I i I i I i I i I i 

-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 

[M/H] 

(b) 

Fig. 10. — Because our Ha index is normalized to the nearby stellar continuum, variability in 
continuum flux as a function of stellar mass or T e ff may affect the interpretation of Figure [9l Here, 
we transform our 7hq values to [Lu a / L^ i\ to remove the influence of the continuum, a. shows 
the time-averaged [Lna/^boi] as a function of stellar mass for our sample. The red line shows our 
linear least squares fit to the data. b. again gives the residuals to the fit shown in a, plotted as a 
function of stellar metallicity. The red line shows our linear least-squares fit to the data. 
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Fig. 11. — Activity levels of M stars has previously been shown to vary with spectral subtype. In 
this Figure, we have grouped our time-averaged [^Ha/^boi] values as a function of spectral subtype 
for our sample. The error bars indicate the RMS scatter in each bin, and the numbers above each 
point give the number of stars in the bin. 
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Fig. 12. — We investigate whether there is any indication that the more metal-poor stars in our sample are older than the metal-rich stars. 
Here, we plot galactocentric velocity dispersion in UVW coordinates as a function of stellar metallicity. Metal-poor stars tend to have higher 
velocity dispersions in all three coordinates, suggesting they are older, potentially explaining why our metal-poor targets exhibit lower mean 
/ho values. 
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Periodic Signals 




Star 


Period (days) 


Amplitude (Iua) 


FAP 


GJ 270 


2687 ± 16 


0.00158 ± 0.0003 


< 10" 4 


GJ 476 


1066 ± 200 


0.00193 ± 0.0005 


< io- 4 


GJ 581 


1633 ± 93 


0.00172 ± 0.0003 


0.12 


GJ 708 


296 ± 1 


0.00127 ± 0.0002 


0.03 


GJ 730 (3-day fit) 


3.06 ±0.01 


0.00128 ± 0.0002 


< 10" 4 


GJ 730 (994-day fit) 


994 ± 40 


0.00134 ± 0.0002 


0.10 


GJ 552 


> 2000 




< 10~ 4 






Linear Trends 








r 


P(r) 


GJ 16 


5.7±i - 2.3 ±0 .5 x W~ 6 t 


-0.868 


0.0012 


GJ 521 


5.4±i - 2.2 ±0 .4 x 10" 6 t 


-0.783 


6.1 x 10" 5 


GJ 96 


2.6 ±0 .5 - l-0±o.2 x 10- 6 t 


-0.720 


5.3 x 10" 6 


GJ 3023 


-8.9± 2 - 3.6 ±0 .7 x W~H 


+0.877 


4.3 x 10" 4 






Quadratic Trends 






lHa(t) 


^poly 


^(^poly) 


GJ 3801 


50000 ± i 6 ooo - 0.04 ±0 .oii + 8± 3 x 10~ 9 i 2 9.31 


0.04 


GJ 611.3 


9000±20oo - 0.007 ±0 .oo2t + 1.5±o. 4 x lO" 9 * 2 12.4 


0.01 


GJ 630 




Not Fit 





Table 1: Parameters of our fits to the In a time series discussed in Section 3. For the periodic fits, 
the FAP listed is taken from the results of our bootstrap resampling analyses. With the linear fits, 
we include the Pearson correlation coefficient r and the resulting probability of no correlation P(r). 
Similarly, the quadratic trends include the results of our F-test, as detailed in Section 3.2.2. Our 
fits to In a (t) are given as a function of barycentric Julian date. 
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Table 2: Ha and Ca I indices shown in Figures H] and [6] 



BJD - 2450000 Ha index Ca I index 

GJ 270 



2619 


74495602 





06040 


± 





00048 





02550 


± 





00034 


2619 


76246122 





06081 


± 





00042 





02548 


± 





00033 


2624 


73217882 





06469 


± 





00054 





02602 


± 





00035 


2635 


93473395 





06218 


± 





00045 





02539 


± 





00033 


2642 


69473208 





05925 


± 





00040 





02588 


± 





00030 


2645 


67800006 





05997 


± 





00039 





02554 


± 





00030 


2649 


66792308 





06027 


± 





00045 





02565 


± 





00035 


2714 


71477040 





06000 


± 





00043 





02510 


± 





00032 


2719 


71740232 





06027 


± 





00041 





02472 


± 





00033 


2756 


61978882 





06236 


± 





00038 





02491 


± 





00032 


2974 


01609612 





06258 


± 





00039 





02576 


± 





00033 


2980 


76850420 





06090 


± 





00040 





02582 


± 





00032 


2987 


75676078 





06166 


± 





00041 





02574 


± 





00034 


2996 


71864146 





06192 


± 





00055 





02574 


± 





00037 


2996 


73387092 





06252 


± 





00039 





02567 


± 





00033 


3003 


71574529 





06248 


± 





00046 





02572 


± 





00033 


3047 


59314927 





06312 


± 





00036 





02536 


± 





00037 


3056 


79170774 





06154 


± 





00066 





02565 


± 





00044 


3070 


74665118 





06062 


± 





00045 





02549 


± 





00034 


3072 


73438156 





06342 


± 





00045 





02500 


± 





00033 


3073 


73829100 





06283 


± 





00045 





02500 


± 





00032 


3096 


66830179 





06060 


± 





00044 





02452 


± 





00034 


3112 


62306948 





06111 


± 





00041 





02495 


± 





00033 


3336 


02956980 





06340 


± 





00045 





02581 


± 





00038 


3358 


75899075 





06146 


± 





00039 





02560 


± 





00038 


3392 


89441676 





06321 


± 





00063 





02631 


± 





00033 


3423 


80129696 





06172 


± 





00043 





02553 


± 





00039 


3704 


00760396 





06272 


± 





00042 





02532 


± 





00034 


3707 


76840342 





06187 


± 





00047 





02571 


± 





00033 


3739 


70652391 





06326 


± 





00050 





02487 


± 





00033 


4193 


67268488 





06071 


± 





00040 





02474 


± 





00033 


4575 


63333239 





06010 


± 





00038 





02528 


± 





00032 


4739 


94479093 





05928 


± 





00042 





02565 


± 





00031 


4811 


76596343 





06041 


± 





00039 





02538 


± 





00033 


5104 


94743407 





05895 


± 





00046 





02568 


± 





00034 


5198 


70136009 





05840 


± 





00035 





02562 


± 





00029 


5241 


58241178 





06112 


± 





00036 





02510 


± 





00036 


5548 


96506060 





06173 


± 





00040 





02536 


± 





00032 


5604 


60284049 





06322 


± 





00037 





02506 


± 





00032 


5847 


91548311 





05990 


± 





00041 





02551 


± 





00030 


5859 


87415031 





06323 


± 





00044 





02537 


± 





00034 


5874 


85664609 





06102 


± 





00042 





02542 


± 





00030 


5875 


84165490 





06175 


± 





00040 





02533 


± 





00031 


5887 


82421737 





05974 


± 





00042 





02535 


± 





00029 


5892 


80041024 





06377 


± 





00047 





02524 


± 





00036 


5895 


78595998 





06432 


± 





00043 





02516 


± 





00032 
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Tabled cont'd. 

BJD - 2450000 Ha index 



5903.76162361 
5916.73727825 
5926.71009501 
5927.70885023 
5953.62135339 
5957.62588797 
5971.59748963 
6031.64221681 
6047.60139792 



0.06174 ± 0.00041 
0.06048 ± 0.00048 
0.06366 ± 0.00037 
0.06478 ± 0.00034 
0.06386 ± 0.00041 
0.06255 ± 0.00037 
0.06217 ± 0.00035 
0.06277 ± 0.00042 
0.06266 ± 0.00038 



Ca I index 
0.02544 ± 0.00031 
0.02547 ± 0.00029 
0.02520 ± 0.00029 
0.02527 ± 0.00030 
0.02491 ± 0.00031 
0.02508 ± 0.00028 
0.02501 ± 0.00040 
0.02412 ± 0.00033 
0.02424 ± 0.00031 



GJ 476 



4809 


00935640 





07405 


± 





00148 





02434 


± 





00088 


4813 


00078938 





07512 


± 





00135 





02449 


± 





00089 


4835 


93825372 





07111 


± 





00138 





02404 


± 





00091 


4840 


93266452 





07162 


± 





00135 





02433 


± 





00089 


4842 


93199363 





06938 


± 





00131 





02414 


± 





00090 


4843 


91050638 





06961 


± 





00134 





02453 


± 





00093 


4845 


91645282 





07082 


± 





00141 





02417 


± 





00092 


5010 


63249904 





07061 


± 





00139 





02457 


± 





00114 


5020 


62065338 





06966 


± 





00134 





02466 


± 





00113 


5183 


98313987 





06975 


± 





00133 





02459 


± 





00093 


5183 


99443947 





06797 


± 





00135 





02460 


± 





00089 


5199 


94793445 


o 


06724 




o 


00136 


o 


02424 




o 


000Q9 


5199 


95853867 





06665 


± 





00142 





02429 


± 





00090 


5207 


92782339 





06894 


± 





00137 





02445 


± 





00097 


5207 


93911867 





06882 


± 





00147 





02419 


± 





00095 


5209 


91363282 





06928 


± 





00148 





02431 


± 





00096 


5240 


99126123 





06805 


± 





00141 





02432 


± 





00090 


5245 


82151704 





06744 


± 





00146 





02442 


± 





00094 


5245 


83282300 





06777 


± 





00139 





02442 


± 





00093 


5253 


95868114 





06881 


± 





00145 





02438 


± 





00094 


5253 


97154526 





06723 


± 





00139 





02445 


± 





00092 


5291 


85653799 





06815 


± 





00140 





02413 


± 





00091 


5291 


86784578 





06702 


± 





00135 





02422 


± 





00092 


5309 


65430856 





06721 


± 





00138 





02424 


± 





00095 


5309 


66572957 





06873 


± 





00135 





02433 


± 





00095 


5337 


73268927 





06683 


± 





00134 





02412 


± 





00097 


5702 


73796815 





06971 


± 





00134 





02411 


± 





00095 


5729 


65679328 





06825 


± 





00139 





02431 


± 





00097 


6017 


87319403 





07135 


± 





00138 





02405 


± 





00091 


6085 


70114197 





06942 


± 





00139 





02401 


± 





00096 


6102 


64742851 





06854 


± 





00139 





02446 


± 





00088 



GJ 581 



2421.74447247 0.07949 ± 0.00156 0.02279 ± 0.00093 
2425.73640218 0.07934 ± 0.00157 0.02230 ± 0.00092 
2443.70418288 0.07827 ± 0.00166 0.02244 ± 0.00099 
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BJD - 2450000 Ha index Ca I index 



2444 


68379899 





07844 


± 





00159 





02212 


± 





00099 


2803 


70286496 





07807 


± 





00159 





02263 


± 





00097 


2806 


69497578 





07957 


± 





00171 





02240 


± 





00095 


2807 


70023704 





08157 


± 





00182 





02258 


± 





00099 


2809 


69168253 





08017 


± 





00176 





02231 


± 





00095 


2822 


66406479 





08091 


± 





00196 





02252 


± 





00099 


2825 


64546199 





08025 


± 





00173 





02237 


± 





00100 


2826 


65113681 





08133 


± 





00182 





02240 


± 





00100 


2832 


64663205 





08170 


± 





00160 





02213 


± 





00098 


3069 


97261868 





07987 


± 





00148 





02188 


± 





00099 


3074 


95259041 





08168 


± 





00143 





02187 


± 





00095 


3433 


98226317 





07716 


± 





00150 





02197 


± 





00103 


3442 


95918646 





07877 


± 





00145 





02220 


± 





00100 


3476 


85336355 





07749 


± 





00134 





02147 


± 





00104 


3481 


84858218 





07608 


± 





00141 





02188 


± 





00105 


3486 


84197848 





07557 


± 





00136 





02202 


± 





00106 


3491 


81260996 





07875 


± 





00142 





02195 


± 





00110 


3507 


79346518 





07902 


± 





00137 





02238 


± 





00104 


3842 


86894559 





07539 


± 





00143 





02236 


± 





00096 


3891 


72779286 





07815 


± 





00134 





02247 


± 





00088 


3891 


74278810 





07925 


± 





00135 





02267 


± 





00090 


4164 


97842690 





08244 


± 





00146 





02207 


± 





00100 


4180 


93453894 





08185 


± 





00144 





02215 


± 





00096 


4217 


83656080 





08090 


± 





00147 





02160 


± 





00109 


4220 


83911388 





07797 


± 





00145 





02198 


± 





00118 


4224 


82363852 





08062 


± 





00145 





02205 


± 





00107 


4550 


92235669 





08065 


± 





00138 





02207 


± 





00097 


4551 


91403140 





07913 


± 





00137 





02218 


± 





00098 


4570 


87339939 





08059 


± 





00138 





02195 


± 





00103 


4575 


84797092 





08029 


± 





00139 





02211 


± 





00101 


4577 


84773598 





08057 


± 





00142 





02217 


± 





00100 


4582 


82691079 





08098 


± 





00144 





02199 


± 





00102 


4609 


77239519 





07976 


± 





00147 





02202 


± 





00090 


4877 


02218442 





07934 


± 





00141 





02190 


± 





00101 


5737 


66382144 





07916 


± 





00173 





02199 


± 





00099 


5740 


66331690 





07898 


± 





00184 





02209 


± 





00100 


5748 


64872936 





07847 


± 





00176 





02204 


± 





00100 


5750 


64396108 





07811 


± 





00187 





02223 


± 





00102 



GJ 708 



2421.76346943 
2421.77979618 
2424.78077980 
2424.95998937 
2441.92344871 
2442.71502056 
2443.90548522 
2816.69540677 



0.05248 ± 0.00052 
0.05280 ± 0.00046 
0.05331 ± 0.00039 
0.05307 ± 0.00037 
0.05121 ± 0.00049 
0.05327 ± 0.00042 
0.05166 ± 0.00041 
0.05171 ± 0.00070 



0.02632 ± 0.00077 
0.02590 ± 0.00073 
0.02596 ± 0.00055 
0.02614 ± 0.00052 
0.02603 ± 0.00040 
0.02597 ± 0.00040 
0.02590 ± 0.00035 
0.02626 ± 0.00066 
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Table® cont'd. 



BJD - 2450000 Ha index Ca I index 



2830 


85081407 





05290 


± 





00071 





02611 


± 





00029 


2832 


86129953 





05139 


± 





00064 





02611 


± 





00025 


3069 


98786632 





05281 


± 





00030 





02518 


± 





00030 


3071 


98819887 





05512 


± 





00037 





02508 


± 





00025 


3073 


96759840 





05579 


± 





00034 





02512 


± 





00027 


3074 


96878585 





05577 


± 





00036 





02419 


± 





00028 


3102 


91202937 





05289 


± 





00036 





02519 


± 





00026 


3157 


95295985 





05190 


± 





00031 





02601 


± 





00036 


3163 


94220809 





05200 


± 





00040 





02633 


± 





00051 


3165 


92941835 





05388 


± 





00049 





02630 


± 





00058 


3172 


91264854 





05093 


± 





00047 





02632 


± 





00046 


3183 


68902746 





05229 


± 





00075 





02633 


± 





00052 


3185 


68858568 





05281 


± 





00088 





02628 


± 





00055 


3188 


87176498 





05318 


± 





00078 





02637 


± 





00043 


3194 


85931903 





05344 


± 





00075 





02601 


± 





00034 


3202 


85173544 





05194 


± 





00073 





02621 


± 





00023 


3205 


84298424 





05290 


± 





00093 





02628 


± 





00023 


3442 


97424666 





05325 


± 





00042 





02575 


± 





00028 


3446 


96351257 





05303 


± 





00031 





02574 


± 





00026 


3451 


94395604 





05277 


± 





00035 





02594 


± 





00025 


3453 


94293104 





05198 


± 





00034 





02572 


± 





00027 


3460 


92021514 





05262 


± 





00031 





02569 


± 





00027 


3474 


87600025 





05147 


± 





00042 





02498 


± 





00031 


3481 


86415042 





05225 


± 





00033 





02550 


± 





00029 


3484 


86326323 





05132 


± 





00039 





02492 


± 





00051 


3487 


86418360 





05134 


± 





00029 





02547 


± 





00032 


3491 


84926346 





05168 


± 





00030 





02547 


± 





00034 


3582 


80419835 





05270 


± 





00071 





02613 


± 





00025 


3590 


76075956 





05143 


± 





00089 





02595 


± 





00025 


3605 


72212062 





05099 


± 





00107 





02614 


± 





00032 


3608 


73157162 





05147 


± 





00088 





02594 


± 





00033 


3652 


58567419 





05410 


± 





00049 





02551 


± 





00026 


3655 


59108942 





05341 


± 





00058 





02567 


± 





00026 


3663 


56322439 





05401 


± 





00060 





02540 


± 





00024 


3668 


57311897 





05268 


± 





00070 





02575 


± 





00034 


3788 


01559027 





05204 


± 





00072 





02547 


± 





00058 


3802 


99870980 





05112 


± 





00062 





02530 


± 





00045 


3806 


97523726 





05163 


± 





00071 





02541 


± 





00040 


3809 


95817637 





05174 


± 





00063 





02522 


± 





00041 


3815 


94969305 





05052 


± 





00052 





02598 


± 





00037 


3898 


72199851 





05059 


± 





00038 





02674 


± 





00053 


3900 


91993971 





05184 


± 





00038 





02633 


± 





00039 


3921 


66287372 





05473 


± 





00123 





02572 


± 





00055 


3939 


82730454 





05080 


± 





00085 





02596 


± 





00023 


3954 


77105354 





05172 


± 





00068 





02558 


± 





00025 


3989 


68370535 





04951 


± 





00080 





02606 


± 





00029 


4161 


99612571 





05024 


± 





00035 





02576 


± 





00030 


4180 


94946715 





05154 


± 





00032 





02522 


± 





00027 


4192 


91495195 





05235 


± 





00031 





02542 


± 





00025 
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Tabled cont'd. 


BJD 


- 2450000 


H 


a index 






Ca I index 




4215 


84574815 





05192 


± 





00033 





02565 


± 





00029 


4230 


81320097 





05202 


± 





00035 





02580 


± 





00043 


4243 


77943634 





05355 


± 





00035 





02583 


± 





00040 


4318 


77353674 





05160 


± 





00103 





02599 


± 





00023 


4340 


72126095 





05073 


± 





00087 





02602 


± 





00025 


4358 


67459053 





05203 


± 





00056 





02601 


± 





00023 


4544 


94380190 





05186 


± 





00031 





02554 


± 





00025 


4574 


87470590 





05325 


± 





00032 





02568 


± 





00025 


4604 


79589362 





05310 


± 





00039 





02582 


± 





00050 


4617 


75601950 





05132 


± 





00037 





02599 


± 





00058 


4632 


71943325 





05107 


± 





00047 





02600 


± 





00052 


4637 


71503991 





05230 


± 





00053 





02603 


± 





00048 


4640 


90649759 





05096 


± 





00059 





02600 


± 





00049 


4908 


95539173 





05186 


± 





00030 





02541 


± 





00025 


5054 


78046302 





05175 


± 





00105 





02602 


± 





00027 


4640 


90649759 





05216 


± 





00042 





02569 


± 





00027 


5801 


72358549 





05069 


± 





00057 





02568 


± 





00025 


5816 


66339932 





05174 


± 





00071 





02560 


± 





00027 


5848 


58886667 





05065 


± 





00042 





02555 


± 





00025 


6002 


95662503 





05173 


± 





00037 





02447 


± 





00028 


6071 


78666575 





05322 


± 





00034 





02561 


± 





00047 


6073 


77503727 





05326 


± 





00044 





02557 


± 





00066 


6082 


95147847 





05295 


± 





00036 





02575 


± 





00046 


GJ 730 


2806 


77210976 





06383 


± 





00065 





02493 


± 





00077 


2807 


78411331 





06422 


± 





00059 





02478 


± 





00076 


2823 


73656281 





06327 


± 





00066 





02452 


± 





00051 


2824 


72252734 





06199 


± 





00077 





02449 


± 





00049 


2824 


73816894 





06226 


± 





00061 





02416 


± 





00046 


2824 


85921550 





06360 


± 





00070 





02458 


± 





00047 


3157 


83409461 





06689 


± 





00062 





02434 


± 





00056 


3182 


87988530 





06254 


± 





00061 





02471 


± 





00059 


3195 


84876306 





06318 


± 





00077 





02473 


± 





00056 


3205 


81091498 





06252 


± 





00098 





02486 


± 





00051 


3215 


78964708 





06323 


± 





00116 





02481 


± 





00041 


3601 


71619604 





06781 


± 





00115 





02440 


± 





00048 


3605 


70795136 





06377 


± 





00126 





02440 


± 





00045 


3609 


71594164 





06570 


± 





00122 





02481 


± 





00048 


3609 


72618630 





06452 


± 





00121 





02454 


± 





00049 


3612 


69830942 





06382 


± 





00096 





02431 


± 





00048 


3615 


70319690 





06379 


± 





00101 





02443 


± 





00049 


3625 


65935946 





06616 


± 





00106 





02430 


± 





00046 


3628 


65472372 





06643 


± 





00089 





02444 


± 





00045 


3899 


77908026 





06165 


± 





00058 





02511 


± 





00069 


3905 


77441403 





06106 


± 





00057 





02482 


± 





00066 


3956 


74984765 





06233 


± 





00121 





02443 


± 





00044 


3966 


72803865 





06378 


± 





00103 





02420 


± 





00045 
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Tabled cont'd. 



BJD - 2450000 Ha index Ca I index 



3997 


64703601 





06266 


± 





00049 





02399 


± 





00045 


4205 


94360136 





06234 


± 





00055 





02444 


± 





00047 


4207 


93884490 





06354 


± 





00055 





02425 


± 





00046 


4249 


95357766 





06339 


± 





00068 





02386 


± 





00063 


4251 


81349095 





06321 


± 





00061 





02403 


± 





00062 


4255 


93527667 





06705 


± 





00067 





02421 


± 





00071 


4329 


61300676 





06384 


± 





00126 





02457 


± 





00045 


4565 


95433274 





06565 


± 





00056 





02380 


± 





00047 


4574 


93691937 





06498 


± 





00061 





02433 


± 





00045 


4586 


90741429 





06459 


± 





00056 





02385 


± 





00048 


4601 


86950800 





06633 


± 





00064 





02372 


± 





00047 


4616 


83025685 





06367 


± 





00075 





02358 


± 





00075 


4663 


70331949 





06496 


± 





00108 





02466 


± 





00073 


4684 


63406145 





06336 


± 





00121 





02439 


± 





00043 


4688 


63578915 





06286 


± 





00136 





02439 


± 





00044 


4732 


62582364 





06284 


± 





00062 





02419 


± 





00045 


4735 


62142047 





06266 


± 





00067 





02417 


± 





00047 


5290 


97110702 





06358 


± 





00054 





02374 


± 





00047 


5290 


97894777 





06366 


± 





00057 





02380 


± 





00047 


5293 


95863282 





06402 


± 





00059 





02396 


± 





00047 


5293 


96646211 





06459 


± 





00058 





02388 


± 





00047 


5294 


96195596 





06410 


± 





00057 





02396 


± 





00047 


5294 


96979451 





06258 


± 





00057 





02398 


± 





00047 


5308 


92314266 





06553 


± 





00056 





02361 


± 





00050 


5308 


93096705 





06555 


± 





00060 





02360 


± 





00047 


5324 


88564373 





06497 


± 





00063 





02331 


± 





00049 


5324 


89347973 





06348 


± 





00066 





02321 


± 





00047 


5351 


81733642 





06626 


± 





00069 





02400 


± 





00066 


5351 


82515888 





06463 


± 





00069 





02412 


± 





00065 


5369 


77055715 





06651 


± 





00066 





02445 


± 





00083 


5369 


77838613 





06435 


± 





00067 





02451 


± 





00081 


5372 


75425063 





06479 


± 





00064 





02448 


± 





00064 


5372 


76207992 





06579 


± 





00060 





02446 


± 





00062 


5800 


71138893 





06316 


± 





00083 





02424 


± 





00048 


5802 


70649889 





06358 


± 





00097 





02437 


± 





00049 


5805 


68692039 





06347 


± 





00103 





02423 


± 





00049 


5808 


68091699 





06326 


± 





00091 





02420 


± 





00047 


5812 


68401024 





06275 


± 





00079 





02421 


± 





00047 


5814 


66188152 





06398 


± 





00084 





02426 


± 





00048 


6064 


86563127 





06316 


± 





00059 





02371 


± 





00045 


du i z 


90U410O I 


u 




± 


A 

u 


uuudo 


A 

u 


Uzo4o 


± 


u 


UUUoU 


6081 


94417888 





06305 


± 





00071 





02375 


± 





00073 


6082 


93894380 





06432 


± 





00072 





02403 


± 





00063 


6088 


92823555 





06389 


± 





00066 





02404 


± 





00062 


6089 


92354455 





06566 


± 





00067 





02432 


± 





00066 


6170 


69649448 





06346 


± 





00079 





02394 


± 





00046 


6170 


71493655 





06120 


± 





00077 





02415 


± 





00047 



GJ 552 
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Tableffl cont'd. 

BJD - 2450000 Ha index 



2448.75113066 
2455.71782301 
2455.73274519 
2469.68925359 
2477.66327025 
2479.66601240 
2480.66247793 
3183.72216292 
3215.63097455 
3215.64649163 
3940.65683466 
4139.90394966 
4218.69833653 
4224.87675309 
4251.79630645 
4460.02422784 
4474.99206740 
4545.81077781 
4560.76617405 
4569.75287007 
4575.92073547 
4584.69254115 
4619.80412414 
4835.00770689 
4859.95361018 
4876.89689302 
5413.62589217 
5695.65820816 



4371.71690550 
4401.79042154 
4402.62106200 
4772.77189820 
5138.77685076 
5542.66322352 
5896.70222250 
5915.64475395 
5929.59721433 



4605.81420282 
4605.82945467 
4606.81447579 
4612.80109233 
4628.73797731 






06696 


± 





00101 





06808 


± 





00100 





06898 


± 





00102 





06694 


± 





00102 





06717 


± 





00103 





06582 


± 





00101 





06572 


± 





00096 





06619 


± 





00100 





06604 


± 





00097 





06618 


± 





00096 





06491 


± 





00092 





06882 


± 





00109 





06925 


± 





00125 





07070 


± 





00117 





06957 


± 





00103 





07025 


± 





00109 





06873 


± 





00107 





07000 


± 





00102 





07106 


± 





00113 





07143 


± 





00108 





07180 


± 





00114 





07083 


± 





00113 





07312 


± 





00116 





06969 


± 





00107 





07314 


± 





00107 





07031 


± 





00111 





06956 


± 





00105 





06948 


± 





00120 



GJ 16 



0.06712 ± 0.00100 
0.06955 ± 0.00091 
0.06685 ± 0.00090 
0.06743 ± 0.00092 
0.06518 ± 0.00091 
0.06451 ± 0.00104 
0.06419 ± 0.00090 
0.06469 ± 0.00099 
0.06475 ± 0.00090 



GJ 521 



0.06859 ± 0.00088 
0.06864 ± 0.00087 
0.06897 ± 0.00086 
0.07064 ± 0.00083 
0.07185 ± 0.00087 



Ca I index 






02412 


± 





00077 





02401 


± 





00075 





02364 


± 





00075 





02404 


± 





00073 





02402 


± 





00075 





02407 


± 





00082 





02396 


± 





00078 





02404 


± 





00077 





02402 


± 





00076 





02395 


± 





00075 





02384 


± 





00074 





02335 


± 





00082 





02379 


± 





00071 





02389 


± 





00074 





02392 


± 





00074 





02390 


± 





00074 





02353 


± 





00077 





02387 


± 





00073 





02380 


± 





00070 





02372 


± 





00077 





02380 


± 





00074 





02384 


± 





00075 





02382 


± 





00077 





02325 


± 





00078 





02364 


± 





00076 





02322 


± 





00082 





02378 


± 





00078 





02358 


± 





00077 



0.02289 ± 0.00082 
0.02364 ± 0.00059 
0.02361 ± 0.00058 
0.02354 ± 0.00057 
0.02351 ± 0.00057 
0.02352 ± 0.00062 
0.02349 ± 0.00064 
0.02355 ± 0.00063 
0.02345 ± 0.00066 



0.02507 ± 0.00075 
0.02524 ± 0.00076 
0.02517 ± 0.00075 
0.02492 ± 0.00101 
0.02451 ± 0.00080 
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BJD 


- 2450000 


H 


a index 






Ca I index 




4813 


02061638 





06894 


± 





00092 





02527 


± 





00065 


4840 


95145462 





06748 


± 





00090 





02523 


± 





00064 


5392 


66208754 





06731 


± 





00081 





02440 


± 





00108 


5684 


88482109 





06586 


± 





00091 





02504 


± 





00067 


5703 


80173796 





06838 


± 





00086 





02490 


± 





00073 


5922 


97859508 





06666 


± 





00094 





02527 


± 





00064 


5947 


91120096 





06510 


± 





00094 





02513 


± 





00064 


6010 


73812761 





06602 


± 





00097 





02516 


± 





00060 


6022 


71712893 





06801 


± 





00099 





02506 


± 





00056 


6078 


77366322 





06744 


± 





00091 





02501 


± 





00068 


6089 


72880678 





06583 


± 





00086 





02449 


± 





00075 


6098 


71812098 





06584 


± 





00088 





02498 


± 





00069 


6110 


68392028 





06612 


± 





00081 





02434 


± 





00090 


GJ 96 


2928 


71737766 





06136 


± 





00044 





02442 


± 





00058 


2928 


72907718 





06212 


± 





00044 





02416 


± 





00056 


2940 


68691338 





06511 


± 





00043 





02392 


± 





00041 


2942 


90722459 





06165 


± 





00048 





02412 


± 





00043 


2944 


67763094 





06087 


± 





00047 





02391 


± 





00045 


2947 


67704335 





06334 


± 





00047 





02347 


± 





00045 


2949 


89161468 





06151 


± 





00052 





02378 


± 





00043 


2951 


63202344 





06196 


± 





00056 





02395 


± 





00042 


2952 


63215148 





06172 


± 





00054 





02371 


± 





00044 


2954 


64677573 





06186 


± 





00057 





02368 


± 





00043 


2958 


84818141 





06215 


± 





00056 





02341 


± 





00041 


2958 


86048483 





06204 


± 





00057 





02325 


± 





00040 


2958 


87197529 





06220 


± 





00058 





02332 


± 





00039 


2963 


64218924 





06372 


± 





00052 





02432 


± 





00039 


2963 


82361268 





06467 


± 





00052 





02446 


± 





00039 


2965 


84855790 





06606 


± 





00054 





02444 


± 





00042 


3295 


68562983 





06174 


± 





00048 





02462 


± 





00057 


3313 


86451512 





06030 


± 





00049 





02475 


± 





00045 


4447 


77277951 





05977 


± 





00052 





02414 


± 





00046 


4727 


99644880 





06175 


± 





00053 





02448 


± 





00048 


5114 


97547929 





05988 


± 





00052 





02446 


± 





00048 


5178 


55049329 





06042 


± 





00050 





02402 


± 





00045 


5546 


56412660 





05914 


± 





00052 





02398 


± 





00045 


5552 


76540903 





06087 


± 





00054 





02354 


± 





00044 


5576 


69626674 





05989 


± 





00062 





02379 


± 





00052 


5863 


90279976 





05977 


± 





00052 





02377 


± 





00043 


5873 


64370409 





05980 


± 





00056 





02353 


± 





00042 


5883 


83733067 





05904 


± 





00058 





02384 


± 





00041 


5910 


53810524 





05935 


± 





00052 





02344 


± 





00045 


GJ 3023 


4341.75211215 


0.06828 ± 0.00089 


0.02267 ± 0.00063 
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Tableffi cont'd. 

BJD - 2450000 Ha index 



4341.76757184 
4366.67935366 
4369.68014645 
4370.68601039 
4482.58914037 
4698.78616303 
5073.75275416 
5567.63575564 
6181.94542826 



0.06663 ± 0.00075 
0.06647 ± 0.00122 
0.06739 ±0.00115 
0.06873 ± 0.00120 
0.06458 ±0.00119 
0.06817 ± 0.00083 
0.06879 ± 0.00099 
0.07211 ± 0.00096 
0.07357 ± 0.00098 



Ca I index 
0.02290 ± 0.00065 
0.02262 ± 0.00105 
0.02347 ± 0.00102 
0.02360 ± 0.00103 
0.02376 ± 0.00070 
0.02295 ± 0.00070 
0.02318 ± 0.00077 
0.02341 ± 0.00072 
0.02179 ± 0.00072 



GJ 3801 



4476 


93940118 





08335 


± 





00187 





02245 


± 





00123 


4476 


95469620 





08605 


± 





00180 





02247 


± 





00122 


4485 


90972080 





07897 


± 





00174 





02321 


± 





00116 


4507 


85448032 





08612 


± 





00172 





02265 


± 





00124 


4518 


82118420 





08597 


± 





00169 





02304 


± 





00120 


4523 


81256612 





08300 


± 





00181 





02303 


± 





00112 


4580 


87665907 





08451 


± 





00181 





02268 


± 





00116 


4628 


76518832 





08241 


± 





00172 





02312 


± 





00120 


4842 


95049497 





08123 


± 





00183 





02288 


± 





00113 


4845 


93432657 





08290 


± 





00172 





02241 


± 





00119 


4869 


86597289 





08198 


± 





00171 





02291 


± 





00125 


5396 


64145117 





07997 


± 





00186 





02263 


± 





00117 


5396 


65275859 





07949 


± 





00188 





02275 


± 





00121 


5685 


64790599 





08241 


± 





00182 





02246 


± 





00119 


6015 


73038865 





08678 


± 





00186 





02261 


± 





00116 



GJ 611.3 



2798 


66234075 





06477 


± 





00062 





02432 


± 





00046 


2798 


68495852 





06285 


± 





00069 





02471 


± 





00047 


2801 


82913296 





06169 


± 





00067 





02474 


± 





00050 


2802 


82265109 





06080 


± 





00073 





02488 


± 





00046 


2805 


81202651 





06222 


± 





00064 





02493 


± 





00050 


2806 


80481737 





06203 


± 





00072 





02463 


± 





00057 


2824 


75534886 





06335 


± 





00075 





02374 


± 





00057 


3136 


92020795 





06245 


± 





00058 





02387 


± 





00060 


3192 


74203561 





06391 


± 





00066 





02359 


± 





00057 


3906 


79408624 





05920 


± 





00054 





02459 


± 





00049 


4575 


95962878 





06157 


± 





00061 





02470 


± 





00052 


4595 


91849395 





06108 


± 





00059 





02395 


± 





00053 


4613 


70348818 





06106 


± 





00078 





02400 


± 





00051 


4616 


85256418 





06125 


± 





00079 





02396 


± 





00048 


4922 


84993631 





06206 


± 





00063 





02466 


± 





00043 


5737 


63188445 





06528 


± 





00061 





02380 


± 





00052 



GJ 630 
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Tableffi cont'd. 



BJD - 2450000 Ha index Ca I index 



2769 


72843715 





05594 


± 





00036 





02634 


± 





00026 


2769 


74429402 





05620 


± 





00033 





02616 


± 





00026 


2770 


74390336 





05811 


± 





00042 





02668 


± 





00026 


2803 


87105450 





05698 


± 





00037 





02666 


± 





00052 


2805 


85276790 





05858 


± 





00057 





02656 


± 





00057 


3138 


95335192 





05525 


± 





00031 





02636 


± 





00024 


3157 


89220154 





05464 


± 





00024 





02662 


± 





00028 


3185 


82205461 





05698 


± 





00045 





02726 


± 





00074 


3202 


78860526 





05915 


± 





00040 





02721 


± 





00057 


3215 


74414667 





05634 


± 





00043 





02740 


± 





00068 


3605 


66947853 





05682 


± 





00051 





02727 


± 





00069 


3608 


65651004 





05565 


± 





00046 





02689 


± 





00065 


3628 


61207839 





05870 


± 





00040 





02723 


± 





00067 


3629 


60629054 





05644 


± 





00035 





02680 


± 





00045 


3640 


58140678 





05911 


± 





00051 





02680 


± 





00059 


3755 


02925396 





05722 


± 





00050 





02683 


± 





00030 


3766 


01195160 





05542 


± 





00047 





02669 


± 





00038 


3770 


00652326 





05820 


± 





00037 





02594 


± 





00028 


3780 


96926864 





05305 


± 





00044 





02783 


± 





00031 


3809 


89071786 





05566 


± 





00053 





02620 


± 





00035 


3816 


88025027 





05509 


± 





00052 





02603 


± 





00033 


3843 


79011280 





05414 


± 





00043 





02664 


± 





00028 


3867 


71826947 





05526 


± 





00024 





02604 


± 





00022 


3872 


93219238 





05430 


± 





00027 





02630 


± 





00024 


3875 


92137530 





05391 


± 





00041 





02585 


± 





00029 


3876 


93499510 





05339 


± 





00026 





02621 


± 





00021 


3879 


93254381 





05458 


± 





00047 





02554 


± 





00026 


3889 


89922516 





05477 


± 





00035 





02639 


± 





00035 


4134 


99091051 





05628 


± 





00051 





02672 


± 





00021 


4136 


98099813 





05487 


± 





00028 





02663 


± 





00019 


4173 


87775570 





05696 


± 





00034 





02635 


± 





00020 


4221 


75954594 





05562 


± 





00027 





02614 


± 





00021 


4251 


67206540 





05652 


± 





00030 





02616 


± 





00030 


4256 


89485013 





05499 


± 





00038 





02586 


± 





00044 


4576 


79143792 





05518 


± 





00051 





02677 


± 





00031 


4588 


75535565 





05501 


± 





00025 





02640 


± 





00022 


4605 


92940501 





05463 


± 





00024 





02610 


± 





00023 


4610 


91897503 





05535 


± 





00032 





02629 


± 





00025 


4617 


89417931 





05460 


± 





00035 





02595 


± 





00024 


4872 


97577035 





05473 


± 





00038 





02665 


± 





00023 


5737 


83092890 





05321 


± 





00033 





02654 


± 





00038 


5737 


84606950 





05358 


± 





00030 





02639 


± 





00043 


5800 


67495127 





05369 


± 





00034 





02704 


± 





00045 


5801 


66454511 





05356 


± 





00033 





02703 


± 





00039 


5806 


64375181 





05345 


± 





00034 





02698 


± 





00051 


5808 


64447093 





05383 


± 





00038 





02687 


± 





00058 


5968 


98110307 





05299 


± 





00028 





02609 


± 





00018 


6014 


83479136 





05560 


± 





00034 





02617 


± 





00022 


6018 


83054636 





05473 


± 





00032 





02636 


± 





00018 
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Tabled cont'd. 


BJD - 2450000 


Ha index 




Ca I index 


6038.78501409 


0.05634 


± 


0.00026 


0.02632 


± 


0.00020 


6063.93996078 


0.05388 


± 


0.00023 


0.02608 


± 


0.00022 


6078.90154931 


0.05573 


± 


0.00027 


0.02613 


± 


0.00024 


6096.85995871 


0.05549 


± 


0.00043 


0.02636 


± 


0.00066 


6145.71936074 


0.05657 


± 


0.00042 


0.02669 


± 


0.00069 


6177.63694756 


0.05735 


± 


0.00051 


0.02687 


± 


0.00080 



